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ABSTRACT 


Hansen 1 s planetary theory represents the position of the disturbed planet as a deviation in 
time and space from the position of a fictitious planet, whose motion is Keplerian relative to 
pseudo-time. This work contains an investigation about the nature of the basic function w of 
Hansen’s theory, which serves to determine the deviation mentioned. The standard way to obtain 
the differential equation for w is based on the application of the method of variation of the astro- 
nomical constants. However, Hansen’s coordinates of a planet are more intimately connected with 
the perturbations in the position vector than with the perturbations in the elements. For this rea- 
son, following von Zeipel, dW/dt is obtained here as the projection of the perturbations in the 
acceleration vector on a variable vector n, which is selected in such a way that dW/dt is totally 
integrable. 

The form of w thus obtained contains an arbitrary constant and an arbitrary function of the 
osculating areal velocity. If dW/dt is obtained in terms of a disturbing force, then it can clearly 
be seen that the presence of the arbitrary elements in w suggests the separate introduction of a 
vector S and of a scalar h 0 /h, bypassing the formation of w. In fact, w represents a fusion of 
three independent series into one, although no actual gain is expected from such a fusion in terms 
of programming or computing time. 

The ”barred” Hansen function w, however, is the most essential part of the theory. It is 
intimately connected with the determination of the perturbations in Hansen’s coordinates and with 
the formation of the integrating operator, w is retained in this exposition. Furthermore, the 
geometric characteristics of Hansen’s theory favor the expansion of the perturbations into trigo- 
nometric series with the disturbed mean anomalies of planets as arguments. Such a form of ex- 
pansion eliminates the need to develop w and the disturbing force in powers of the perturbations 
of the mean anomalies. It contracts the series and speeds up the convergence. We economize, 
especially, in the expansions of the odd negative powers of the mutual planetary distances. 

The motion of the disturbed planet is related to an inertial frame of reference situated in the 
orbital plane of the fictitious planet. The standard use of the rotating ideal frame of reference 
instead of an inertial frame does not simplify the actual numerical procedure. Hansen himself was 
compelled to introduce the perturbations in the ’ ’third coordinate.” In the final instance this is 
equivalent to the use of a fixed frame of reference of the type mentioned. 

The integration operator is expanded into a series, whose form favors the application of the 
method of iteration. At each step of the iterative process the perturbations are obtained by solv- 
ing a linear partial differential equation with constant coefficients. 
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NOTATION 


P - the disturbed planet, 
m - the mass of P, 

M - the mass of the Sun, 
f - the gravitational constant, 

/a 2 = f(M+m), 

Pj - the auxiliary planet moving in a fixed plane in accordance with Kepler's laws, 

M, a o> e o> * * * etc. the standard elliptic elements of P 1? 

P' - the disturbing planet, 
s - the heliocentric position vector of P, 

7 - the projection of s on the orbital plane of P, 
r ~ \ 7 I , 

r - the heliocentric position vector of Pj , 

7 z I f | , 

7 - the true anomaly of P x , 

E - the eccentric anomaly of Pj , 
u - the distance of P from the orbital plane of P x , 

P - the unit vector directed toward the perihelion of P x , 

R - the unit vector normal to the orbital plane of P, , 

Q □ R xP, 

l + v - the ratio r/r, 
t - time, 

z - the pseudo-time (the disturbed time). It appears instead of t in Kepler's equation for P x , 
d7 

v D dz~ " the velocity of r relative to z, 

8z = z - t - the perturbations of time, 

n 0 8z - the perturbations in the mean anomaly, 

The notations for P' corresponding to the above are designated by the "primed” symbols, 

J - the mutual inclination of the orbital planes of two auxiliary planets, 

V' - the gradient operator relative to 7' , 

X - 7 ' - 7 - the mutual distance between two auxiliary planets, 

87 = vf + uR - Hansen perturbations of 7, 

7 + 87 = 7, 

87' - v'r'+u'R' - Hansen perturbations of r', 
r ' = r ' + 8 r ' , 

D ' - V exp (87 ■ V' ), 

D " = V ' exp (S'' ■ V ' ) , 

I - the operator of projection on the orbital plane of P ls 

l/h - the areal velocity of the projection of P on the orbital plane of P x , 

l/h 0 - the areal velocity of P 1? 

X = h 0 /h, 

W and w - Hansen's functions. 
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INVESTIGATIONS IN HANSEN’S PLANETARY THEORY 


by 

Peter Musen 

Goddard Space Flight Center 


INTRODUCTION 

This work presents an investigation concerning the foundations of Hansen 1 s planetary theory 
(1857), and a system of formulas which can be used in actual computations of the first and higher 
order general perturbations in Hansen's coordinates. 

As the first approximation to the motion of the planet P we take, following the suggestion by 
Oppolzer (1883) and Andoyer (1926), the motion of a fictitious auxiliary planet P x moving in a con- 
stant ellipse in a fixed plane (xy) in accordance with Kepler T s laws. The position of P is determined 
by its deviation from the position of P 2 in time and space. 

More precisely: at the moment t the projection of the position vector of P on the (xy)-plane 
will have the same direction as the position vector of P 1 at the moment z. Designating by s the 
heliocentric position vector of P, by r the heliocentric position vector of the projection of P on (xy), 
by u the distance of P from (xy), by r the heliocentric position vector of and finally by R the 
unit vector normal to (xy), we can write the basic relations of the theory in the form: 

s(t) = r(t) + u(t) R , 

r(t) = (l+^)7(z), (1) 

The fictitious planet should be so chosen that Hansen's perturbations u, and Sz in the semi- 
major axis, in the "third coordinate", and in time, respectively, will be small— of the order of the 
perturbations. The introduction of the pseudo-time z and of the "perturbations of time” Sz is the 
main characteristic of Hansen's theory which was considered peculiar during his lifetime and 
which, surprisingly, is sometimes considered difficult to accept even now. Another important 
characteristic of Hansen's theory is that the perturbations v and Sz in the orbit plane of P x , are 
determined by one single function W. The mean anomalies of the fictitious planets corresponding 
to the disturbed and to the disturbing planets are the basic arguments in expansions of the per- 
turbations into periodic series. The role of the mean anomaly of P 2 is twofold: It enters into the 
expansion of the perturbations of Hansen's coordinates and it is also the argument in the expansion 
of the elliptic coordinates of P,. In order to separate the perturbative effects from the elliptic 
motion of Pj Hansen uses a special notation for each mean anomaly. Then he forms a function W and 
its derivative dW/dt in which these anomalies are kept separated. The derivative dW/dt is formed 
in terms of the components of the disturbing force. The W-function is obtained by integrating 
dW/dt. After the integration, the W-function is obtained by removing the distinction between the two 
types of mean anomalies. 
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The operator for converting W into W is Hansen’s ’’bar-operator," and its use represents the 
third main characteristic of Hansen’s theory. Normally the differential equation governing the 
variation of W is deduced on the basis of the variation of elements. Thus, the theory of Hansen can 
be considered as an ingenious transfer of the perturbations in the elements to the perturbations 
in the coordinates. However, an approach without resort to the variation of elements and based 
solely on the equations of motion is also possible. Such an approach was suggested by von Zeipel 
(1902). He used the Lagrangian form of the equations of motion in polar coordinates with the eccen- 
tric anomaly as the independent variable. 

Looking closely at dW/dt, one can recognize easily that it is the projection of the perturbations 
in the acceleration of P on a variable vector selected so that dW/dt is totally integrable and w is 
a linear function of dv/dt. W contains an arbitrary vector c which can be taken to be a function of 
the ’’elliptic” mean anomaly. Hansen’s solution corresponds to the special case when c is normal 
to the position vector of the auxiliary planet. We follow here in von Zeipel’ s footsteps with some 
modifications. The direct use of Equation (1) leads immediately to the expression for the pertur- 
bations of the acceleration vector of P and to the general forms of dW/dt and w which facilitate 
the kinematical conclusions. 

Looking at these general forms of w and dW/dt we recognize that w, rather than w, represents 
the main feature of the method. The W-function has a direct kinematical meaning, whereas w repre- 
sents merely an artificial device to combine these series into one. At the present time the ex- 
pansion of the perturbations into series is done on electronic machines. Neither the length of the 
series, nor the computing time can be reduced by employing W instead of the three series. For 
this reason we propose to compute the three parts of w separately. After the computation is com- 
pleted, we form the W-function. 

We deviate from Hansen in the method of expanding the disturbing force and in the method of 
integration. In Hansen’s work the expansion of the disturbing function and of the disturbing force 
is done in powers of v , v', n 0 Sz, n 0 ' 8z', u and u'. The expansion in powers of v, v ' and u, u' is as- 
sociated with the expansion of the disturbing force in odd negative powers of the mutual distance of 
the planets, while the expansion in powers of n 0 Sz and n 0 ' Sz' is associated with the expansion of the 
odd negative powers of l as functions of time. The perturbations n 0 Sz and n 0 ' Sz' normally are the 
largest ones. Their determination is associated with the possible appearance in the process of 
integration of the squares of the small divisors caused by the commensurability of the mean motions 

Thus, the convergence in powers of n 0 Sz and n 0 ' Sz' is a relatively slow process as compared 
to the convergence in v, v', u, u'. To speed up convergence we discard the development in powers 
of n 0 Szand n 0 ' Sz I Then the angular arguments in the periodic series representing the perturba- 
tions or their derivatives will become the linear functions of the perturbed mean anomalies M and 
M'. The integration of such series is then reduced either to solving a partial differential equation 
if we use the process of iteration, or to solving a chain of linear partial differential equations, all 
of the same type, if we resort to the standard asymptotic expansion of perturbations in powers of 
masses. 
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We will not follow in Hansen’s and von Zeipel’s footsteps and do not employ the eccentric an- 
omaly as the independent variable. However, only a slight modification is required in the method 
we present if it is desired to switch to either the eccentric or the true anomaly as the basic variable. 

The computational scheme is arranged so that the process of iteration is applicable and thus 
the progra mm ing is homogeneous. If preferred, the computational scheme for the standard asym- 
ptotic expansion in powers of masses can be easily deduced. 

THE DIFFERENTIAL EQUATIONS OF MOTION AND THE EXPANSION 
OF THE DISTURBING FORCE 

The motion of the planet P t is Keplerian and its position vector 7(z) satisfies the differential 
equation 


d 2 r _ _ fj?_ r 

dz 2 " 7 3 


( 2 ) 


The equation of motion of P can be written in the form (Musen, 1965) 


A 2 

dt 


377 (7 +uR) 


(r 2 +u 2 ) 


3/2 


( r + uR) + fm' 



(3) 


where 


l = l?l . (4) 

l * 7' - 7 , (5) 


D'' = V' exp(sl • V') , 

D' = V' exp ( S? • V') , (6) 

S-t — v' ~f‘ ~ v * + u' R' — uR = 8 r' -Sr , (7) 


Sr - v r + uR , 

Sr" 1 = v' T' + u' R' , (8) 
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and V' is the gradient operator relative to r'. From the author’s previous (1965) work, we take the 
expansions: 


D " X 


+ {is 11 '* 1 - l S ^ 1 + 


15 

V 


+ — - 1 * Si 8l + l Si • hi 
l 5 l b 


1(1 ■ Si) 2 


+ 1+ w - s?) 3 - £ a • S^) 2 S? 


|| + • Si 


+ 1(1 • 8 ?) (5-e • 8-t) + — 84 • 54 84 ^ + • * * ; (9) 


and 


D' 


/ 3 — 

- l 

-\ r 


'•87'- — 

8r') + [- 

6 


3 n ■ 

+ 7^ 

r "F ; • 8 r ' Sr' 

+ t= 7 j r 8 

r 3 

45 

r 


“ T 77 

(?' * 87' ) 2 7 

' +7' (7' 


•' • S7')j +^j (87' • 87') Sr 


• (10) 


In the process of expanding (9) and (10) and the components of the disturbing force in terms of the 
disturbed mean anomalies, we shall use the following auxiliary quantities: 


B = B' = 


a ■ "»■ (f) ■ 

A ■ 

• r ' = a Q a Q ' ( acc # + /3cs / + a' sc' +/3' ss ' ) 

F = r' ■ R = a Q ' (yc ' + y' s') , 

F' = r * R' - a Q (kc + k' s) , 

A' = R * 7 x 7' = a Q a 0 ' (a' cc ' -asc' +/3' cs' -/3ss') , 


(11) 

( 12 ) 

(13) 

(14) 

(15) 
(13’) 
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( 14 ') 

( 11 ") 

(13") 

(15") 


c 

-u° 

II 

:os f , 

s 

r 

sin f , 


c ' 

r ' 

cos f ' , 

s ' 

r ' ' 
_ / 

sin 7', 



a o 



a o 



a 

= P • 

P' , 

a ' 

= Q • 

P' , 


0 

® P • 

Q' , 

£' 

= Q • 

Q' , 


y 

= R • 

P' , 

y' 

= Q' 

■ R , 


K 

= P • 

R' , 

K ' 

- Q • 

R' . 



G 

- r 

• l = 

< 

1 

PQ 

* 

(16) 


G' 

- 7' 

II 

A' - B' . 

(17) 


The expansion of (11) - (17) in terms of M and M' is performed in the standard manner, either using 
the classical analytical formulas or by employing harmonic analysis. 

The explicitly written portions of (9) and (10) permit the development of the general perturba- 
tions up to the fourth order relative to the disturbing masses. We have from (7): 

2 • M = G' v' - Gv - F' u' - Fu , (18) 
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S-C • S-t - A' v ' 2 + Av 2 - 2Bv ' v - 2FV 'u - 2F ' vu' + (u 2 + u ' 2 - 2uu ' cos j) 


7' ■ 87' = A' v' , 


87' • 87' = A' v' 2 + u' 2 . 


(19) 

( 20 ) 
( 21 ) 


Substituting these expressions into (9) and (10) we deduce the expansions: 


D " i + D' 


Kr + K' 7' + HR + H' R' 


( 22 ) 


K 


/ * K .C n .X / i/i u k u /k ; 

4 i j . J ,k,k 


(23) 


C* - / K'.( n ). » / u k u' k ; 

/ . j , j ,k, k » 


(24) 


H - / * H / v* v‘ * u k u ,k 

/ , J > J * k , k 


(25) 


H' - ^ * j/j v* * u k u' k , 

^ , j.j *k,k 


(26) 


j , j ' , k, k ' ^ 0, j + j ' + k + k ' 


n = 0, 1, 2, 


In the following table we give the values of the coefficients K and H which are different from 
zero for n = 0, 1, 2. 


K (°) = - — 

^0000 £3 ’ 


K ' <°> 
^oooo 


_1 1_ 

/ £>3 7' 3 ’ 


K 

000 


JL _ 3G 

t 3 t 5 9 


K ' (1 > 
1 ooo 


3G 

l 3 9 
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K <*> 
0100 


3G ' 

I s 


K ' O) 
*0100 


1 3G / 2 

-e 3 ~ -e s + ?'3 


K O) 
ooi o 


3F 


K ' (!) 
^ooz 0 


3F 
< s • 


K 

0001 


3F' 

-e s ’ 


K ' (!) = + — 

0001 + 


H ofto 


-e 3 ’ 


H 0 0 ( 01 


I _ JL_ 

4 3 7/3 * 


K (2) = _ 15. r 2 6G , 3A 

20^0 ^7 G - ^ 


K (2) 
0200 


— — q ' 2 + — — — 

• 


K (2) 
0020 


15 XT' 2 3 

F 2 + — 

V t* • 


k o ( o 2 } 2 


- — p' 2 + ± 

V is ’ 


K < 2 ) = + — nr' 4 . 6G " 6 B 

1100 47 GG + Jf ~ j? , 


K ( 2 ) 
010 1 


30 , 

+ ^ G F ' . 


K (ni 0 = 


+ 30 6 F 

^ fg <ts . 


K (2) 
1001 


30 . 12F' 

■ — GF 

V V 


K (2) 
*1 0 1 0 


30 rp .6F 
't 7 -e* - 


K o ( o 2 L 


30 , 6 cos J 

V is 


7 



K 7 (2) 
^0200 


6 r , 3 A' 

4 * G 4 s 

K # < 2 > 
^2000 

■ 

3 A 

~ 4 s ' 

K 7 < 2 ) 
^0002 

■ 

4 s 7 ' 5 

K ' (2) 

"“01 1 0 

30 ^ 12 F 

" v G F+ -e 5 ’ 

K 7 < 2 > 
^2 001 

30 

+ — F' 

4 s 

K ' < 2 ) 
*1010 

30 


K ' < 2 > 
R ooi 1 


6 

+ v cos J ■ 

H (2) 
n 01 1 0 

6G 7 

+ t 5 ’ 


H < 2 > 
n i 01 0 

= -H 7 ( 2 ) = 

1001 

6G 
4 s ’ 

H < 2 > 
n ooi 1 

= -H ' ( 2 > = 

n 0 0 0 2 

6F' 

4 s ’ 

n 002 0 

= -H '< 2 > = 

001 1 

6F 

4 s ’ 


H '(2) = - + -i!- . 

n 0101 45 T J '3 

This table permits one to develop the perturbations up to the third order. This accuracy is 
sufficient for the majority of planets in our solar system. The inclusion of the perturbative ef- 
fects of the fourth and higher orders, if necessary, can be easily accomplished by expanding the 
additional terms in (9) and (10). 

The choice of the orbital planes of the auxiliary planets can be made in such a way that the 
orders of the magnitudes of u and u ' numerically will be higher than the orders of magnitudes of 
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V, v' , n 0 Sz, n 0 ' Sz'. Thus, the terms in F which contain the factors m'vu, m'v U ', m' u' u , m' v' u 
m' u 2 , m' uu', m' u' 2 can usually be omitted. 

From (3) we have 


d 2 r 
dt 2 



(27) 


d 2 u 
dt 2 " 


u 


M - 

* 2 + u 2 ) 3/2 


+ z 


(28) 


where we set 


F - M 2 


(r 2 + u 2 ) 


Z - R 


2 \ 3/2 


r + 


fm'I • ^-D" j + D' ^ 


- D' j t D 1 r 
^ r ' 


(29) 

(30) 


I is the operator of projection on the (xy)-plane. It can be represented in several ways, either as 

I ” h Q Rx(rv-vr) 


h Q ( r v - v r ) x R , 


(31) 


or as 


I - ttj (rf-RxrrxR) , 


(32) 


or as 


I - PP + QQ , (33) 

according to the kind of decomposition of F we prefer. The representation (31) of I is preferable 
if we want to separate the terms with small divisors. The first terms in (29) and (30) can be 
written as 


(r 2 +u 2 ) 3/2 


3 9 u 2 

T ^ =s C 1 + 


(34) 
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and 


2 

jJL U 


(r 2 +u 2 ) 


3/2 


fJL 2 U 3 7 u 2 

T + 2^ + 


The last quantity normally is negligible. 

Taking (11) - (15") into account, we deduce the following expressions for Z and for the com- 
ponents of F: 


Z = fm' (K' F + H +H' cos J) 


(35) 


7 ■ F = fm' (KA + K' B + H' F') + ^ p 2 (l+^)~ 4 + 


R x 7 • F = fm' (K' A' + H' ®' ) , 


=• 3 u 2 

v * F = fm' (KT+ K' F' +H' $' ) + *2 /x 2 — (1 + ^)' 4 F + 


(36) 

(37) 

(38) 


P • F 


3 u 2 

fm' [Ka 0 c +K' a 0 ' (ac' + /3s ' ) + H' k] + ^ ^ 2 a Q — ( 1 + )' 4 c + 


(36’) 


•3 m ^ 

Q-F = fm' [Ka 0 s + K' a 0 ' (a' c' + /T s') + H' #c '] + \ (x 2 a Q ~ ’ ( 1 + v)‘ 4 s + (37 f ) 


The system of formulas (35) - (37 T ) ? because of its symmetry, provides us with a scheme which is 
convenient for the use of the process of iteration as well as for programming. 


PERTURBATIONS IN THE ORBITAL PLANE OF THE AUXILIARY PLANET 

Differentiating (1) twice with respect to time and taking (2) and (27) into account, we have 


v - vr + (l + i^)vz 


(39) 


T.. M 2 ( 1 + v) z 2 l ^ V d r _ . o o 

73 — j r + rrv dt [c 1 + v ) z ] 




+ F 


( 40 ) 
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Introducing Hansen's notations for the sectorial velocities, 


R • r x 7 - h 


R • r x v 


(41) 

(42) 


we deduce from (39) the classical formula 


(1 +v) 2 


(43) 


where we set 


TT = x • 


(44) 


Taking (43) into account we can rewrite (40) as 


x2 " (! + vS ) tL 

(1 + v) 3 7 3 


r + 


1 + v 


(45) 


We now determine a vector N so that the scalar product of N with the left side of (45) represents a 
totally integrable function. 

For this purpose a convenient form of the decomposition of N is 


N - pvxR + qRxr 


(46) 


where p and q are functions of X, v, z. From (42) we have 


p - h 0 r • N 


q - h Q v ■ N 


(47) 

(48) 


We determine a function 


W = W(X, v, v, z) 


(49) 
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such that, in accordance with (45) - (48), we have 


dW 

dt 


i? . iw .. 

d* x + di , v 


'aw . aw A 
v + a7 z ) 



X 2 - (1 + v) y?_ q 1 
(1 + v) 3 . T 3 ] + h o 1 + v 


( 50 ) 


In this respect we follow in von Zeipel’s footsteps. However, since we are using time, and not the 
eccentric anomaly, as the independent variable, the straight forward application of (2) and (27) 
makes the exposition simpler. In addition, the kinematical meaning of the results becomes more 
direct. 

From the last equation we deduce: 


aW _ _1_ q 

ax h 0 i + v » 


(51) 


aw = _p 

dv h o 


(52) 


aw . aw 

dv v + dz 


X 

(1+v) 2 


(i +*) ^ x 2 . jSp_ m 
(i+vy t 3 h o 


(53) 


From (52) we conclude that W must have the form 


= ^ £ + b (X , L*. Z ) 


Substituting (54) into (51) and (53) we obtain 


ap 

a^7 = 0 3 


*p _ 0 

ax " 0 


ab / 1 x \ /x 2 p 
az ~ \x ” i +v) -3 ’ 


_ab x dp 

dv ~ (1 + v) 2 dz 


ab _ q 

ax ~ i + v * 


(54) 


(55) 

(56) 

(57) 
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Thus p is a function of z only and from (55) and (56), we deduce 


d 2 p _ _ m 2 p 

dz 2 F 3 


( 58 ) 


The general solution of (58) is 


p - h 0 c • r , 


(59) 


where c is an arbitrary vector lying in the orbital plane of Pj. From (59), we also deduce 


dp _ - = 

dz " h 0 c • v 


From (55), (56) and (60), we have 


(60) 


b 


( JL 

1 0 \ 1 + J 



V + h 0 f(X) 


(61) 


where f(X) is an arbitrary function. 
From (57) we have 


q = h o (l + c • F + h 0 (l + v)f' (X) 


(62) 


Substituting (61) and (62) into (46) and (54) we obtain 


N - h Q c * rv X R + h Q 


7 1 + zA _ =* 

(1 + J C • V+(l +!/) f 


R x r 


(63) 


and the general form of the W -function is 


W = v! • c + 




(64) 


Taking (31) into account, we obtain 


_ / ho I- A „ 

N - c • 1 1 + Rx rl + h Q f R x ,r 


( 65 ) 
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From (45) and (50), we deduce 


dW - - 

dF = N • F . 


Substituting (64) and (65) into the last equation, and because c and f are arbitrary, we obtain 


dX - ^ - 

dF = h 0 R • r x F , 


S' r + (iru - x) v ] = F + X 2 V R • r 


x F 


(66) 

(67) 


The equation (66) is classical. The equation (67) is new. It is convenient to bypass the use of w 
and introduce the vector S defined by the differential equation 


dS 

dt 


h 0 R x 


(f ?R- ?xf) . 


( 68 ) 


Of course, this formula can be connected with the variation of astronomical constants. 
We have from (67) 

1 7 + (rTU ~ x) v = ixRh " 1 . 


(69) 


Taking (42) into consideration we deduce from the last equation, after some easy vectorial 
tr ansf or mations , 


= + v 


(70) 


X i _ — ~ 

1 + v ~ X " " r ' b 


(71) 


We can obtain v from (70) by integration. Then a constant of integration will appear which de- 
pends upon the constant in X. Perhaps it is preferable to determine v from (71) which does not in- 
volve any additional integration. We can re-write (71) in the form 


V 


(1 +v) 



which is convenient for the application of the process of iteration. 


(72) 
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At this step we introduce Hansen's W-funetion in order to determine n 0 Sz. We have 


w 


1 - X + 


2X 

1 + v * 


The quantity X differs from unity by the order of perturbations and we can set 


X = 1 + 8X . 


( 73 ) 


(74) 


Using (71) we obtain 


W = - (l + SX - 2? • S , (75) 

and the perturbations in the mean anomaly and z can be obtained by using Hill's (1881) formulas: 


(76) 

(77) 

PERTURBATIONS IN THE THIRD COORDINATE 

Taking into account the identity 

d 2 u = Td 2 u du dt d /dtYj /dz\ 2 , 

dt 2 " " \d^~ " dt dz dt \ dz /J\ dt / ’ 

we can rewrite (28) in the form: 


dt 


+ v* 


dz 

dt 


1 + W 


d 2 u u 2 u dt 

~ + c di 


(78) 


where we set 


1 dz 

1 dt 


du 

d /dt\ 

—3 dt 

(r 2 + u 2 ) 3/2 dz _ 

u + Z dl 

+ ar 

dt (dz/ 


(79) 
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Expanding the first term in powers of u and using (71), (43) and (77), we obtain 


C 


1 - V 2 _ 2 1 U - 

1 + W ^ - /x 1 +^ _ 3 r 


S + 


du d 
dt dt 



( 80 ) 


From (43) we have 


d /dt\ _ 2( 1 + v ) d^ ( 1 +v) 2 dX 

dt \dz ) X dt ^2 dt 


From (73) we obtain with the needed accuracy 


1 + V 

X 


1 - \ (W + SX) 


(81) 


and substituting this value into (81) we have 

Tt (ai) = (2-W-8X) f t - (1-W-SX) f- ■ 


Thus (80) becomes 


*^Z- M ’(l-„) •— ? • S + 
1 + W ^ K T 3 


__ dv — dX 

(2-W-8X) ^ " (1-W-SX) dF 


du 3 , u£ 

dt + 2 M 2 ?5 


(82) 


For the majority of planets the perturbations of the second order in u will suffice, providing that 
the orbital plane of P' is chosen in a proper way. Then we can put 


_ ( dv dX\ du \j? u =t - 

c = (l-W)Z + (2 dt - dFjdF-^ r ■ s • 


(83) 


The integration of (78) gives 


U 



C r p sin (<£ 


f)dt 


(84) 


At this point it was convenient to use the Hansen and Hill device: the auxiliary mean anomaly (M) 
in the expressions of 

■§- cos 4 > and sin<£ 
a o o 
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I 


is considered as a temporary constant. After the integration is completed (M) must be replaced by 
M. In computing the terms of the second order in (83) the elliptic values of the coordinates can be 
used. 

THE PROCESS OF INTEGRATION 

The integration of the differential equations of the problem is accomplished by the method of 
successive approximation using the expansion of the integrating operator into series. The operator 


d 

D = dt 

can be decomposed into a sum of two operators 

D = T 0 + T 


(85) 


where 


_d_ , _d_ 

T o ” n o dM + n o 


(85 ? ) 


and 


T 


d ( n o^f) d_ d ( n o' s f^) 
dt c?M + dt 


d 

dM' 



* 1/12 d 

u 2 dM + ! - „'2 dU > 


(85”) 


If only the perturbative effects of the first and second orders are needed then we set v 2 - v ' 2 - o 
in (85*) and the operator T takes a simpler form. Putting (85) in the form 


D = Mlo + T- 1 !-) , 


where l 0 is the identity operator, we obtain the expansion of the integrating operator. 


Q = D - 1 = (i 0 +t- 1 t)-»t 0 -i 


T ” 1 _ T-l TT“1 + T“1 TT“1 TT“1 - ... 

1 0 1 0 T1 o 1 A o AA o 


( 86 ) 
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The explicitly written terms in this expansion suffice for computing of the perturbations up to the 
third order. The central part of the integration procedure is the application of the integrating 
operator T 0 '" 1 . The basic formulas to be used are either 


To" 1 ( zS z' s ‘) 


1 — S+ 1 's _ •** „s + 2 ' s ' - 1 . ~ „ s +3 's -2 

S+1 Z Z (S + 1)(S+2) Z Z + + Z 


_jL( s ' - 1) 


(s +1) (s + 2) (s +3) 


s ' + 1 


, s 7 1 s +1 __ 


s(s - 1) 


(s' +1) (s' + 2) 


(s' + 1) (s' + 2) (s' +3) 


? S“2 is +3 _ 


and 


V 1 {z* «'■' exp(jM + j'M')} = y _^}± k + J_V (2 S z 's' } 

0 l u 1 J n 0 + j n 0 4^ (i n 0 + j n o ) k V dz ) K 1 

or their modif ications . In practical application of these formulas, we usually do not go beyond 

m = 2. 

After v u, n 0 Sz are obtained, the time is determined from the equation 


n 0 t - n Q z + c Q - n 0 Sz - M - n 0 Sz . 


If the t is given then the corresponding value of z is determined from the last equation by means of 
successive approximations. The disturbed position vector s is determined using the standard 
equation 


= ( 1 +V) 


a Q P (cos E~ e 0 ) + a 0 - eT 0 2 Qsin E + uR 


CONCLUSION 

The geometrical characteristics of Hansen's planetary theory favor the expansion of perturba- 
tions in terms of the disturbed mean anomalies of the auxiliary planets. This approach eliminates 
the need to expand in powers of n 0 Sz, n 0 ' SzJ * • . and speeds up the convergence of series repre- 
senting Hansen's coordinates. We economize especially on the expansions of odd negative powers * 

of the mutual distances. Thus, we are not compelled to form the derivatives of the first and 
higher orders of the disturbing function relative to the pseudo-time. Hansen's theory is connected 
more intimately with perturbations in the rectangular coordinates than with perturbations of the 
elements. For this reason we suggest here the form of the theory which does not appeal to the 
method of variation of astronomical constants. We prefer instead to make direct use of the dif- 
ferential equations of motion in rectangular coordinates relative to an inertial frame of reference. 
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We did this because the classical use of the moving ideal system of coordinates evidently does not 
simplify the actual numerical procedure. 

Hansen himself was compelled to introduce the perturbations of the "third coordinate. " To 
obtain them he used a method which combines the direct kinematical considerations with the method 
of variation of constants. Furthermore, we do not make use of Hansen’s W-function because of the 
way the operations with periodic series are performed on electronic machines. However, we con- 
sider employing W as the most essential part of Hansen’s theory. The W-function is intimately 
connected with the determination of the perturbations in the mean anomaly as well as with the 
formation of the integrating operators. For these reasons we have retained the W-function in the 
present exposition. 
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